Creating two-dimensional bright solitons in dipolar Bose-Einstein condensates 
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We propose a realistic experimental setup for creating quasi-two-dimensional (2D) bright solitons 
in dipolar Bose-Einstein condensates (BECs), the existence of which was proposed in Phys. Rev. 
Lett. 100, 090406 (2008). A challenging feature of the expected solitons is their strong inher- 
ent anisotropy, due to the necessary in-plane orientation of the local moments in the dipolar gas. 
This may be the first chance of making multidimensional matter-wave solitons, as well as solitons 
featuring the anistropy due to their intrinsic dynamics. Our analysis is based on the extended Gross- 
Pitaevskii equation, which includes three-body losses and noise in the scattering length, induced 
by fluctuations of currents inducing the necessary magnetic fields, which are factors crucial to the 
adequate description of experimental conditions. By means of systematic 3D simulations, we find a 
ramping scenario for the change of the scattering length and trap frequencies which results in the 
creation of robust solitons, that readily withstand the concomitant excitation of the condensate. 

PACS numbers: 03.75.-b, 03.75.Lm, 05.45.Yv, 67.85.-d 



Bose-Einstein condensation (BEC) of magnetic atoms, 
or molecules carrying electric moments, has attracted a 
great deal of attention, starting from the pioneering the- 
oretical analyses of such condensates [l[ , which was fol- 
lowed by the creation of BEC in the ultracold gas of 52 Cr 
atoms, using magnetic Q and all-optical Q techniques. 
Very recently, a condensate of 164 Dy atoms, with a much 
larger magnetic moment, was created too jjj. The mak- 
ing of magnetic- and electric-dipolar BEC may also be 
expected, respectively, in erbium [H[ and in molecular 
gases Q . An updated review of this rapidly progressing 
field was given in Ref. Q- 

One of the promising possibilities, which has been, 
thus far, investigated only theoretically, is the creation of 
bright solitons in dipolar condensates. This can be easily 
predicted in various one-dimensional (ID) settings, us- 
ing, in particular, periodic potentials induced by optical 
lattices as the stabilizing factor 0] . In the limit of a very 
deep lattice, the BEC wave function becomes nearly dis- 
crete, which makes it possible to predict stable 2D dipolar 
solitons and vortices [l(|, as well as 3D solitons frij . 
in the (quasi-) discrete form. 

As concerns solitons, the most challenging issue is the 
creation of such stable matter-wave modes in quasi-2D 
(pancake-shaped) condensates, as well as their counter- 
parts in the form of spatiotcmporal solitons in optics. In 
spite of intensive theoretical discussions (l2j . no experi- 
mental results in this area have been reported thus far, 
the most essential obstacle being the inherent instability 
of 2D solitons to the collapse, driven by the self-focusing 
cubic nonlinearity (very recently, it was proposed to cre- 
ate stable bright solitons by means of self-defocusing non- 
linearity, which is possible in a rather exotic situation 
with the strength of the nonlinearity growing atr^oo 
faster than r D , where D is the spatial dimension [HI). 



The dipolar condensates suggest new challenges and pos- 
sibilities for achieving this fundamental purpose, by mak- 
ing use of the competition between local interactions and 
long-range dipole-dipole interactions (DDIs). A quasi-2D 
isotropic configuration implies that the local moments are 
polarized perpendicular to the pancake's plane, in which 
case the DDI is repulsive. In that case, the creation of 
(bright) 2D solitons may be possible if the sign of the DDI 
is effectively reversed by means of a rapidly oscillating 
magnetic field 14J; in the same setting, stable isotropic 
solitons with embedded vorticity were predicted too [15| , 
and various 2D localized structures may be stabilized by 
trapping potentials acting in the plane [lfj . The very fact 
of the existence of multidimensional solitons in the dipo- 
lar BEC can be proven in a rigorous mathematical form 
17j |. It is also relevant to mention that stable isotropic 
solitons are possible in optical media with nonlocal (ther- 
mal or orientational) nonlincarities, as was demonstrated 
experimentally [18| . and studied in detail theoretically 
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The most challenging issue, which is unique to the mul- 
tidimensional dipolar media, is the possibility of the cre- 
ation of a new species of anisotropic solitons, based on 
the in-plane polarization of local moments (obviously, the 
anisotropy cannot manifest itself in ID solitons). The 
stationary form of such quasi-2D solitons, and some of 
their dynamical properties, were investigated in Rcfs. 
[20L l2lT |. but the problem of the actual creation of the 
solitons under experimentally relevant conditions has not 
been addressed before. The complex anisotropic struc- 
ture expected in the solitons makes this problem crucially 
important, as the mode may be considered physically ro- 
bust only if it is capable to self-trap from experimentally 
feasible initial configurations. In this work, we show, by 
means of systematic simulations of a realistic model, that 
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anisotropic solitons in dipolar BEC can be created, us- 
ing available experimental techniques. This result may 
have more general implications, suggesting the realizabil- 
ity of various complex modes predicted in multidimen- 
sional quantum gases, such as skyrmions [22T |. 

As mentioned above, there are currently two atomic 
species available for creating magnetic dipolar BEC, viz., 
52 Cr and 164 Dy. Since the latter has been condensed 
very lately Q, our main results refer to chromium, but 
at the end of this work we address dysprosium too. 
The making of the anisotropic solitons is possible in a 
strongly dipolar regime, which may be defined by the 
dimensionless ratio of the DDI and local interaction, 
£dd = ra^iQ^L 2 / il2-n:afi 2 ). Here, m is the atom's mass, 
Ho the free-space permeability, \i the atom's magnetic 
moment, and a the s-wave scattering length. As shown 
in Ref. [20|, the necessary condition for the existence of 
solitons is 6dd > 1- In the case of 52 Cr, the natural 
value of this ratio is 6dd ~ 0.16 [2^]. Therefore, it is 
necessary to tune the scattering length to enhance the 
relative strength of the DDI, which may be accomplished 
by using the Feshbach resonance [3,[2^]. For 52 Cr, we 
thus find that the critical scattering length below which 
the 2D solitons can exist is 15 as- However, this condi- 
tion is not sufficient, as it is only an upper bound, and 
the soliton may actually be found at still smaller val- 
ues of a 2l|. On the other hand, if a is too small, i.e. 
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are the trap frequencies in the respective di- 



£dd exceeds a certain critical value, the condensate will 
collapse. Therefore, solitons only exist within a certain 
range of values of the scattering length. For typical trap 
frequencies and numbers of atoms in the condensate, the 
relevant range for chromium is in a ballpark of 4 as 21 1 . 

To explore conditions for the creation of the anisotropic 
solitons, we consider BEC of N atoms with magnetic 
dipole moment fj,, which are fully polarized by means of 
an external uniform magnetic field. At a sufficiently low 
temperature, the dynamics of the condensate obey the 
extended Gross-Pitaevskii equation (CPE): 



h (t) + font (t) + feioss(i) #(r,t) = ihdt^(r,t), (1) 



with the single-particle Hamiltonian 



ho(t) = ~A + y [c4(i)z 2 + W 2 y 2 + u? z {t)z 2 \ 
the mean-field interaction Hamiltonian 
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rcctions, L3 determines the rate of three-body losses, and 
\£ is the mean-field wave function. The dipoles are polar- 
ized along the z-axis, and it is assumed that the trap fre- 
quencies in x- and z-directions may be time-dependent, 
along with the scattering length of the contact interac- 
tion, a. We aim to consider the following scenario: Start 
with the dipolar BEC at parameters for which the con- 
densate is stable in the extended state filling the entire 
3D trap, and ramp down the scattering length and fre- 
quencies uj x , uj z simultaneously to values at which a sta- 
ble anisotropic soliton is expected. In the course of do- 
ing this, the trap will be switched off completely in the 
x- and z-directions. Obviously, this condensate-steering 
scenario will excite the condensate. We will exploit this 
fact for proving the self-trapping nature of the DDIs, 
since stable oscillations without the external trap are 
only possible if the long-range interactions provide for 
the self-confinement. To make our model realistic, we 
include losses via the imaginary-potential term (j4]) , with 
an empirically known loss coefficient, L3 = 2 x 10~ 40 
m 6 /s [2^, [27j], and random noise in the scattering length, 
which stems from the current in the magnetic coils in- 
ducing the Feshbach resonance. The latter is of special 
significance, because a typical root-mean-square value of 
the noise, ~ 1 — 2 ob , is not much smaller than the range 
of a within which the soliton exists. The noise that we 
used in the simulations is a uniform distribution of all 
frequencies up to 500 Hz. We do not expect higher fre- 
quencies to significantly affect the dynamics of the soli- 
ton because the time scale of these dynamics is ~ 1 s, 
as shown below. We have analyzed the noise signal of 
the Feshbach-inducing current in the experimental setup 
employed at the laboratory in Stuttgart (28|, observing 
close agreement with the uniform distribution adopted in 
our scheme. 

The ramping sequence is shown in Fig. Q] The ini- 
tial condensate contains 20000 atoms, and is prepared 
in the trap with (uj x ,oj y ,oj z ) — 2tt x (72, 426, 72) Hz and 
a = 27 as ■ The scattering length is then ramped down to 
11.5 ob within 70 ms. The trap is first held constant for 
thold = 4 ms and is then ramped to f X)Z = fi = 6.5 Hz 
within 58 ms. Finally, it is switched off within subsequent 
33 ms. We emphasize the necessity for eliminating the 
trap. Even a very weak trapping potential, with / ~ a 
few Hz, affects the behavior of the soliton, making it very 
difficult to distinguish the resulting dynamics from that 
of an ordinary trapped condensate. For the simulations, 
the 3D spatial domain was discrctized into 512 x 128 x 512 
grid points, and marching in time is performed using the 
Split-Operator method. The DDI integral is evaluated 
with the aid of the Fourier convolution formula. Since 
this numerical scheme is very demanding, the code was 
implemented using the parallel computing architecture 
CUDA, and it was run on a modern graphics card, en- 
abling a very high degree of parallelization. 
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FIG. 1. (Color online) Ramping of (a) scattering length a and 
(b) trap frequencies. Only the first 100 ms are shown. In the 
course of the remaining time of the simulations, the scattering 
length is held at a constant value, while the trap has already 
been switched off. In the experiment, the scattering length 
cannot be held at an exact value, but contains some noise, 
which is included into our scheme, as shown by the green line 
in panel (a). 



Results of two such simulation runs are displayed in 
Fig. O The left column shows absorption images of the 
condensate (|^| 2 integrated over the y-direction and nor- 
malized to the maximum value) in the absence of the 
noise in the scattering length, while, in the right col- 
umn, the noise with the root-mean-square value of 0.5 Ob 
was added. In the first two frames, the soliton is being 
shaped. It is observed that few atoms escape the core 
condensate in the x-direction, along which the DDI is 
repulsive. The condensate then grows further in the x- 
and z-dircctions, but eventually inverts its dynamics and 
rcfocuses. Without the noise, this occurs at < w 303 ms. 
Since this is long after the trap has been switched off, the 
self-confining potential needed for this outcome can only 
be provided by the DDIs. Thus, we find a clear evidence 
for the self-trapping of the robust soliton. Remarkably, 
the condensate refocuses as well when the noise is added 
to the scattering length, as is shown in the right column 
of Fig.f2J i.e., the buildup of the soliton is robust against 
this experimentally inevitable noise. 

To gain more insight into the dynamics, we monitored 
the root-mean-square extensions of the condensate in the 
x- and z-dircctions during the self-trapping process, as 
shown in Fig. [3l It is observed that the expansion in 
these directions is roughly synchronous, which implies 
that a breathing mode has been excited. 

When the noise is added to the scattering length, the 
condensate reaches a greater extension, and the oscilla- 
tion frequency decreases. Above a critical amplitude of 
the noise, for which we can give the lower bound of 1 ob , 
we expect the soliton to decay. This behavior can be un- 
derstood in terms of a simple picture: One can replace 
the noise by an effectively larger scattering length, which 
also lowers the binding energy of the system and there- 
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FIG. 2. (Color) The formation of a soliton. Left column: 
The scattering length is ramped without the noise. Right 
column: The random noise with the root-mean-square value 
of Aa — 0.5 Ob was added to the scattering length. Shown 
are absorption images of the condensate density, with \^>\ 2 
integrated along the y-axis and normalized to the maximum 
value. In both cases, the soliton is being shaped in the first 
two frames, and is still expanding at t = 217 ms. Without the 
noise, the condensate self-focuses already at t = 303 ms, while 
this is slowed down when the noise is included. In both cases, 
the refocusing occurs at the evolution stage when the external 
trap has already been switched off completely, hence it is the 
DDIs which are responsible for the self-trapping. This is a 
clear indication of the existence of a stable excited soliton. 
It is important to note that the noise does not destroy the 
soliton, which is crucial for the experiments. The field of the 
view is 232 pm x213 /im. 
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FIG. 3. (Color online) Root-mean-square extensions of the 
soliton in (a) x- and (b) ^-directions for various values of 
the noise. Naturally, the largest extension of the condensate 
increases when the noise is added to the scattering length. At 
the end of the simulations performed without the noise, the 
soliton has already started a second cycle of the oscillations. 
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fore slows down its dynamics. Consistent with Rcf. [2l| . 
the increase of the scattering length, a, also leads to a 
larger extension of the soliton. Above the critical value 
of a, the soliton decays into a freely expanding conden- 
sate. From the experimental point of view, the noise may 
even help to watch the soliton dynamics in situ because 
of the comparably large observable change in the conden- 
sate extensions. 

The small dip in the x-extension at around 100 ms is 
noteworthy, as it is caused by the still present external 
trap with a small frequency, /i, which decelerates the 
expansion. In line with this trend, increasing f\ from 6.5 
Hz to 8 Hz we have observed a more pronounced dip and 
a smaller value of the maximum extension, together with 
an increased oscillation frequency. For smaller holding 
times ihoid, meaning that the ramping sequence for the 
trap begins earlier, we find a shallower dip, leading to 
slower dynamics and larger extensions of the condensate. 

The simulations show that the time scale governing 
the dynamics of the soliton is on the order of 1 s, which 
is relatively large for trapped condensates. The main 
limitation for holding times is imposed by losses caused 
by three-particle collisions, which is modeled in Eq. ([T]) 
by imaginary potential (U|) proportional to the squared 
number density. From Fig. [31 we observe that, in the 
case of the ideal ramp, the soliton is larger in the x- and 
z-directions by a factor of 4— 5 than the initial condensate 
(in the y-direction, the size is approximately constant). 
This results in the density reduced by factor ~ 15 — 25, 
hence the loss rate will rapidly drop when the trap is 
switched off. Indeed, Fig. [4] where the evolution of the 
number of atoms in the condensate during the simula- 
tion is plotted, demonstrates that most atoms arc lost 
within the first 50 ms. Without the noise in the scatter- 
ing length, the atom number then stays nearly constant 
for subsequent 300 — 400 ms, and decreases at the end 
of the simulation, when the soliton has reached its small- 
est extension, i.e., it has returned to the largest density. 
With the noise present, the mean extensions are even 
larger, therefore the losses arc reduced even more. For 
the 52 Cr condensate with the scattering length around 
1 1 q,b , we thus conclude that the atom losses will not 
cause decay of the soliton. 

Finally, for parameters of the 164 Dy condensate we find 
that the threshold scattering length necessary for the ex- 
istence of the anisotropic solitons is 133 as- Since the 
actual scattering length of Dy is estimated to be smaller 
than this value |4| , anisotropic solitons may be created in 
this condensate by means of a simpler procedure, gradu- 
ally switching off the trap only and keeping the scattering 
length constant. Applying a simple Gaussian trial wave 
function to evaluate the mean-field energy of a conden- 
sate which is only trapped in y-direction, as was done 
in Refs. 0, tHJ, we can estimate the binding energies 
of the soliton in the (x, z) plane for various scattering 
lengths. This energy scale determines the frequency of 
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FIG. 4. (Color online) The number of atoms in the conden- 
sate in the course of the formation of the soliton. Most atoms 
are lost during the first 50 ms, when the density is still large 
and the soliton has not yet been formed completely. When 
the soliton starts its first cycle of oscillations, the loss rate 
rapidly drops, and at the end of the simulation, the conden- 
sate still keeps more than 80% of the initial atoms. For the 
noise amplitude of 1 as , the losses are not visible at all after 
the trap has been switched off. 



oscillation of the excited soliton, showing a crucial de- 
pendence on the scattering length. For the confinement 
of 20000 atoms in the y-direction, we choose a trap fre- 
quency of 200 Hz, and conclude that the relevant time 
scale varies from 0.5 s to 2 s for the scattering length 
between 126 cir and 122 Ob- At 128 ob, the time scale 
already reaches a value ~ 10 s. 

In conclusion, we have performed systematic simula- 
tions of the realistic model describing the evolution of 
the dipolar condensate en route to the formation of the 
recently predicted new species of localized multidimen- 
sional solitary modes, in the form of quasi-2D anisotropic 
solitons. The model includes factors which are crucially 
important to the adequate description of the experiment, 
such as the three-body losses and Feshbach-induced noise 
in the scattering length. The scenario for the forma- 
tion of the solitons was elaborated, by gradually switch- 
ing off the trap and decreasing the scattering length. It 
is demonstrated that, in the chromium and dysprosium 
BEC alike, the solitons may be safely formed, surviving 
the inevitable excitation of the condensate. Further chal- 
lenging problems may be the making of soliton clusters, 
and the study of interactions between anisotropic soli- 
tons. A more general implication of the reported results 
is the possibility to elaborate similar dynamical scenar- 
ios for the creation, under realistic conditions, of vari- 
ous complex multidimensional patterns predicted by the 
analysis in quantum gases. 
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